Simulation of current-induced microwave oscillation 
in geometrically confined domain wall. 
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We studied magnetization dynamics of a geometrically confined domain wall under dc current by 
solving simultaneously the Landau-Lifshitz-Gilbert equation and diffusion equation for spin accumu- 
lation. We showed that the oscillation motion of the domain wall is driven by the spin-transfer torque 
and the dc current is converted to the ac voltage signal. The results means that the geometrically 
confined domain wall is applicable as a source of microwave oscillator. 



Achievement of a nano spin-transfer oscillator is one of 
important issues in microwave technology. Many studies 
have succeeded in development of microwave oscillation 
in nanopillers and point contacts [l|, [3, d, 0, 0, 0, 0, H] ■ Re- 
cently possibilities of another kind of nano spin-transfer 
oscillators using a domain wall were proposed^, [lol |. 
In those works, the oscillation generating microwave in 
GMR or TMR devices is not an uniform oscillation of the 
free layer magnetization but an oscillation [Q] or magne- 
tization rotational motion [lo| of the domain wall. 

On the other hand, it is known that the resistance due 
to the domain wall depends on its magnetic structure pdl 
m, [ill • If th*^ oscillation motion of the magnetic struc- 
ture of the domain wall is induced by spin transfer torque 
under dc current, we can obtain the ac voltage signal and 
can use the domain wall as a microwave source. 

In this paper, we numerically calculate the dynamics of 
the geometrically confined domain wall under the dc 
current by solving simultaneously the Landau-Lifshitz- 
Gilbert equation and diffusion equation for spin accumu- 
lation. We show that the oscillation of the magnetization 
vector of the domain wall is induced by the spin-transfer 
torque and we can obtain the ac voltage signal due to 
the domain wall oscillation. We also show that the power 
spectrum of the MR ratio has a sharp peak correspond- 
ing to the oscillation of the domain wall. The frequency 
of the peak is about 1 hundred GHz for the current den- 
sity of about 0.01 mA/nm^ which is reachable value in 
experiments with avoiding breakdown by Joule heating. 

The system, we consider, is shown in Fig. [T] The sys- 
tem consists of two magnetic electrode layers and a non- 
magnetic insulator layer sandwiched by the electrodes. 
The insulator layer contains magnetic bridge which geo- 
metrically confines a domain wall. The system is divided 
into about 1500 hexahcdral finite elements, where spin 
accumulation and demagnetization field are evaluated. 
The size of magnetic metal electrode at the bottom and 
top layers is 8 x 8 x 4 nm. The shape of the bridge 
between two electrodes is a rotated elliptic arch around 
the center axis perpendicular to the layers, the diameters 
of the bottom and center layers of which are set at 6 and 
4 nm, respectively. The size of the system is larger than 
the diffusion length of 2nm of CoFe[la| at room temper- 
atures. Thus we can deal with the electron system in 
diffusive limit, unlike a previous ballistic treatment . 
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FIG. 1: The model of a geometrically confined domain wall 
for a finite element method. 



Then as previously predicted by Simanekp^. the spin 
accumulation mainly induces domain wall resistance. 

The magnetization vectors arc expressed by the clas- 
sical spins on the simple cubic lattice with the lattice 
constant of a = 0.4nm. The Hamiltonian Ti. is given by 

H = - Jdd ^ Si ■ Sj + Jsd ^ Si ■ Srhi 



where represents the relative coordinate of the i-th 
site from the position r, S'i, the classical Heiscnberg spin 
with absolute value of unity and, Srfii^ local spin accu- 
mulation density at the i-th site. Such the semi-classical 
treatment is justified for the system with ferromagnetic 
order because the magnetization is larger than its quan- 
tum fluctuation. The first and second terms in the right 
hand side of Eq. ^ express the exchange interactions be- 
tween local magnetizations at nearest neighbor sites, and 
between local magnetization and spin accumulation at 
each site, respectively. The exchange coupling constant 
between nearest local magnetization denoted by Jjd is 
related to the exchange stiffness constant in continuous 
limit, A, and a lattice constant, a, with Jdd = 2aA. We 
set Jdd at 0.04 eV which is on the order of the transition 
temperature of CoFe, Tc ^ 1.2 x IQi^K. The exchange 
coupling constant Jgd is set at 0.1 eV in accordance with 
Ref . p7j . The third term denotes the dipole-dipole in- 
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teraction energy. In order to deal with the dipole-dipole 
interaction in a complex system shape, we adopt a fi- 
nite element - boundary element (FEM-BEM) hybrid 
method [isj on the spin field on the continuum space, 
5(r), which is defined for each position, r, by the spin 
field extrapolated from the lattice sites to the position. 
The exchange length for the dipole-dipole interaction is 
Zex = \/ Jdd/'^Kd ~ 3 nm, which is comparable to the 
length of the considered confining region of 4nm. Thus 
the thickness of the domain wall depends on both the 
size of the confining region and the exchange length. 

In the adiabatic approximation, the local spin accumu- 
lation density, Srrii , is determined by solving the following 
differential equations , 



^(5m(f) = V {/35(f)/o(r) + A(S{f))5m{r)'^ 

-I- —5m{r) X Sir) -\ 

n T 

A{S{r)) = 2L>o \i - [3'^S{r) ® S{f)v\ , 



(2) 
(3) 



where je, Cq, -Doj P a-nd r denote electronic current, con- 
ductivity, diffusion constant, polarization of resistivity 
and relaxation time due to a spin-orbit interaction re- 
spectively. Equations (HI)-® are solved numerically with 
combining continuous equation for electronic current, 
Vje = 0. Co and (3 are taken to be those for conventional 
ferromagnetsl as Co =70rinm and (3 = 0.65, respectively. 
Dq and r are obtained, respectively, by the Einstein rela- 
tion, Co = 2e^NFDo and A = y/2TDoil-P^) for given 
diffusion length A, density of states at the Fermi level, 
A^F and electron charge e. Here we employ A=12nm and 
Np =7.5nm~^eV~^. For the boundary condition of the 
spin accumulation, we artificially adopt Srn = at the 
top and bottom layers, ignoring any parasitic resistance. 
On the other boundaries, the natural boundary condition 
is employed. In the case, the simulated spin accumula- 
tion distribution is nonuniform and concentrate on the 
contact region. 

Substituting the solution of Eqs. ^ and ^ into 
Eq. ^ , we evaluate the Hamiltonian of Eq. ([T|) . The ef- 
fective magnetic field for Si is given by dTi/dSi and the 
dynamics of Si is determined by the following Landau- 
Lifshitz-Gilbert equation. 
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were 7 and a denote the gyromagnctic ratio and the 
Gilbert damping constant. The equation is numerically 
solved by a quaternion method [l9[. The time step At 
and a are set at 3.4 x 10~^ fs and 0.02, respectively. 
The number of sites, Ng is about 10^. For simplicity, we 
use antiparallcl boundary condition of the spins in order 
to simulate the experimental situation with a 180° do- 
main wall. On the boundary, S's are fixed at (-1,0,0) on 
the top electrode and (-1-1,0,0) on the bottom electrode 
except for the boundary surface between the magnetic 
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FIG. 2: (a) The x, y, and jz-components of the mean mag- 
netization vector, Af, is plotted as a function of time t for 
j = 0.012mA/nm2(b) The power spectra of the magnetiza- 
tion, Px,y,z, are plotted as a function of frequency u). 



and insulator layers. The directions of spins, x, y and z 
are defined as shown in Fig. [TJ 

In the present work, we employ current of 0.012 
mA/nm^. The current is applied in the perpendicular 
direction to the layers. The simulation starts from sta- 
ble Neel wall state obtained by a simulation without dc 
current density. 

In order to understand dynamics of the local magne- 
tization distribution of the domain wall, we calculate 

the magnetization M = {Ylii Sij /Ng for the whole of 

the system. The typical result in time dependence of 
the magnetization is shown in Fig. [HJa). We observe 
steady regular rotation around the x-axis. Such rota- 
tion is known for the moving domain wall in wire sys- 
tems above Walker's threshold current |2(|. The magnetic 
structures with My—0 and M^—O correspond to the Neel 
wall and the Bloch wall, respectively. Thus the rota- 
tion corresponds to an oscillation between the Neel and 
Bloch walls. On the other hand, the displacement of the 
domain wall center is characterized by the value of 
because the magnetization vectors at the boundary sur- 
face arc aligned to be parallel to the x-axis. As shown in 
FigEJa) that the domain wall displacement oscillates as 
a function of time. 

Figure [2jb) shows the power spectra of the magnetiza- 
tion motions, which are defined by 
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(T)Mx,y,,(t-T)dTdi, 



(5) 

where Tq and T are much longer than the periods of the 
oscillations. Each spectrum has a single characteristic 
frequency. The frequency of the displacement oscilla- 
tion is twice of that of the rotational motion. This is 
because the displacement oscillation originates from cou- 
pling with the rotational motion driven by the spin trans- 
fer torque. For the Neel and Bloch walls, respectively, the 
displacement takes minimum and maximum. 

The oscillation of the magnetization vectors induces 
the oscillation of the spin accumulation through Eq. 
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FIG. 3: (a) The magnetoresistance ratio Rmb. is plotted as a 
function of time t for O.OlmA/nm'^. (b) The power spectrum 
of the magnetoresistance ratio Pmr is plotted as a function 
of frequency ui. 



The electric field at position r is given by 



2Do 



i + l3a-S{r) Srh{f)] . (6) 



The voltage drop of the system is obtained by integrat- 
ing the electric field E{r) along z-axis. From Eq. the 
oscillation in the spin accumulation then yields the oscil- 
lation in magnetoresistance ratio, i?MR(i) = R{t) /Ro — 1, 
where Ro and R{t) denote resistance with and without a 
domain wall, respectively. Rmr oscillates as a function 
of time t as shown in Fig. [S^a). The power spectrum of 
-Rmr defined by 
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PuRitu) = ^ JJ^ e-''^'RMRiT)RMR{t - T)dTdt (7) 

is shown in Fig. (S^b). In the accuracy of our calculation, 



Pmr has a characteristic isolated peak. Thus the dc cur- 
rent is converted to the ac voltage signal with a charac- 
teristic frequency. The frequency of the peak agrees with 
that of the denoted by the arrow in Fig[2)^b). Thus the 
MR oscillation is due to periodically change of the spin 
accumulation induced by the displacement oscillation. 

We note that only the resistance due to the spin accu- 
mulation is considered. However in the real system the 
resistance originates from not only the spin accumula- 
tion but also scattering of electrons due to the domain 
wall [Hill. The correction will enhance the magnetore- 
sistance ratio and therefore the power of the microwave 
oscillation as compared with the results of our simulation. 

In conclusion, we studied the magnetization dynamics 
of the geometrically confined domain wall under dc cur- 
rent by solving the Landau-Lifshitz-Gilbert equation and 
diffusion equation for the spin accumulation. From calcu- 
lation result, we propose a scenario of the microwave gen- 
eration as follows: current induced spin-transfer torque 
drives magnetization rotation in the domain wall. Then 
the displacement oscillation of domain wall is induced by 
the coupling with the rotation and drives spin accumu- 
lation oscillation. As a result the microwave oscillation 
in voltage signal appears. We conclude that the geomet- 
rically confined domain wall is applicable as a source of 
microwave generator. 

The authors would like to thank M. Doi, H. Iwasaki, 
M. Ichimura, M. Takagishi, M. Sahashi, M. Sasaki, 
T. Taniguchi, N. Yokoshi and K. Seki for useful dis- 
cussions. The work was supported by The New En- 
ergy and Industrial Technology Development Organiza- 
tion (NEDO). K. M. was supported by a Grant-in- Aid for 
Young Scientists from the Ministry of Education, Science, 
Sports and Culture of Japan. 



[1] J. A. Katine, F. J. Albert, R. A. Buhrman, E. B. Myers, 
and D. C. Ralph, Nature 84, 4212 (2000). 

[2] M. Tsoi, A. G. M. Jansen, J. Bass, W.-C. Chiang, 
V. Tsoi, and R Wyder, Nature 406, 46 (2000). 

[3] S. I. Kiselev, J. C. Sankey, I. N. Krivorotov, N. G. Em- 
ley, R. J. Schoelkopf, R. A. Buhrman, and D. C. Ralph, 
Nature 425, 380 (2003). 

[4] W. H. Rippard, M. R. Pufall, S. Kaka, S. E. Russek, and 
T. J. Silva, Phys. Rev. Lett. 92, 027201 (2004). 

[5] M. Govington, M. A. H. Darwish, Y. Ding, N. J. Goke- 
mejier, and M.Seigler, Phys. Rev. B 69, 184406 (2004). 

[6] I. N. Krivorotov, N. G. Emley, J. G. Sankey, S. I. Kise- 
lev, D. G. Raph, and R. A. Buhrman, Science 307, 228 
(2005). 

[7] S. Kaka, M. R. Pufall, W. H. Rippard, T. J. Silva, S. E. 

Russek, and J. A. Katine, Nature 437, 389 (2005). 
[8] F. B. Mancoff, N. D. Rizzo, B. N. Engel, and S. Tehrani, 

Nature 437, 393 (2005). 
[9] J. He and S. Zhang, Appl. Phys. Lett. 90, 142508 (2007). 



[10] T. Ono and Y. Nakatani, Appl. Phys. Exp. 1, 061301 
(2008). 

[11] P. M. Levy and S. Zhang, Phys. Rev. Lett. 79, 5110 
(1997). 

[12] A. Brataas, G. Tatara, and G. E. W. Bauer, 60, 3406 
(1999). 

[13] E. Simanek, 63, 224412 (2001). 

[14] P. Bruno, Phys. Rev. Lett. 83, 2425 (1999). 

[15] S. M. Moyerman, J. G. Gheckelsky, S. S. Harberger, A. G. 

Tamboh, M. J. Garey, P. D. Sparks, and J. G. Eckert, 97, 

10G513 (2005). 
[16] J. i. Ohe and B. Kramer, 96, 027204 (2006). 
[17] S. Zhang, P. M. Levy, and A. Pert, Phys. Rev. Lett. 88, 

236601 (2002). 
[18] D. R. Fredkin and T. R. Koehler, 26, 415 (1990). 
[19] B. Visscher and X. Feng, 65, 104412 (2002). 
[20] G. Tatara and H. Kohno, Phys. Rev. Lett. 67, 113316 

(2004). 



